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Abstract. The Joint Replenishment Problem (JRP) is a fundamental 
optimization problem in the supply chain management concerned with 
optimizing the flow of goods between a supplier and retailers. In JRP, 
retailers receive demands for a commodity. To satisfy these demands they 
place orders at the warehouse which, in turn, issues aggregate orders 
at the supplier. The objective is to compute a schedule of orders that 
minimizes the sum of the ordering costs and delay costs. 
In this paper, we study the approximability of JRP-D, the variant of 
JRP where there is no delay cost, but the demands needs to be satisfied 
before specified deadlines. The main result of the paper is an e/(e — 
l)-approximation algorithm based on iterative rounding, improving the 
previous bound of 5/3. For the case of equal-length demand periods, we 
improve the ratio to 1.5 and we show that this case remains APX-hard 
even if each retailer has at most four demands. 

1 Introduction 

The Joint Replenishment Problem with Deadlines (JRP-D) is an optimization 
problem in the inventory theory concerned with optimizing the flow of goods 
between suppliers and retailers. JRP-D arises in a distribution setting with one 
supplier, one warehouse and m retailers, cf. Figure 1. The retailers receive de- 
mands for a certain commodity that they need to provide within a specified 
time period. To serve these demands, the retailers place orders at the ware- 
house which, in turn, places orders at the supplier. As any single order comes 
with a cost, it is desired to aggregate the orders at each retailer and at the 
warehouse, while still assuring that the goods reach retailers in time. 

More concretely, the problem specification involves a fixed warehouse order- 
ing cost L and, for each retailer p, his retailer ordering cost bp. A demand is 
represented by a triple (p, r, d), where p £ {1, . . . , m} is the identifier of the re- 
tailer, r G Q is the demand's release time and d G Q is its deadline. The interval 
[r, d] will be referred to as the demand period. At any chosen time t, a subset R 
of retailers may place an order, through the warehouse, at the supplier, which is 
then served immediately and satisfies all the demands (p, r, d) , such that p € R, 
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Fig. 1. An illustration of an instance with four retailers, where the ordering costs are 
represented graphically as distances. The cost of an order is the total weight of the 
subtree connecting the supplier and the involved retailers. 

and t S [r, d] . The cost of this order is then equal to the ordering cost of the 
warehouse plus the ordering costs of respective retailers, i.e., L + X^^efl ^p- It 
convenient to think of this order as consisting of a warehouse order of cost L, 
which is then joined by each retailer p G Rat cost bp. A schedule is a collection of 
orders. A schedule is feasible if it satisfies all the demands. The objective of the 
JRP-D problem is to find — for any given set of demands — a feasible schedule 
that minimizes the total cost. 

Previous Results. The decision variant of JRP-D problem was shown to be 
strongly NP-complete by Becchetti et al. [3] . (They used a different terminology, 
considering an equivalent problem of packet aggregation with deadlines on two- 
level trees.) Later, Nonner and Souza [10] showed that JRP-D is APX-hard, 
even in the special case when each retailer issues only three demands. A 2- 
approximation LP-based algorithm was given by Levi, Roundy and Shmoys [7]. 
The approximation ratio was subsequently improved to 1.8 by Levi et al. [6,8,9], 
using randomized rounding, and then to 5/3 by Nonner and Souza [10]. 

Our Contribution. The main contribution of our paper is to further improve 
the approximation ratio to e/ (e— 1) ~ 1.58. Our algorithm is based on a recently 
developed technique of iterative randomized rounding of an LP-relaxation, thus 
showing a more general applicability of this approach. We also consider a spe- 
cific case in which all demand periods are of equal length. From the practical 
standpoint, this would capture the scenario where the demand period lengths 
are set globally, for example via industry standards. We show that this version 
is as hard as the general case, by proving that JRP-D remains APX-hard with 
this restriction, in fact even if each retailer has only four demands. On the other 
hand, for equal-length demand periods we show a simple 1.5-approximation. 

Related Work. JRP-D is a sub-case of a more general Joint Replenishment 
Problem (JRP). In JRP, instead of having a deadline, each demand is associated 
with a delay cost function that specifies the cost for the delay between the times 
the demand was released and satisfied by an order. JRP is NP-complete, even if 
the delay function is linear [2,10]. JRP is in turn a sub-case of One- Warehouse 
Multi-Retailer (OWMR) problem, where the commodities may be stored at the 
warehouse for a certain cost per time unit. For instance, the 1.8-approximation 
result by Levi et al. [9] holds also for the general setting of OWMR. JRP was 



also studied in the online scenario, where a 3-competitive algorithm was given 
by Buchbinder et al. [4] . 

Another straightforward generalization of JRP involves a tree-like structure 
with the supplier in the root and retailers at the leaves. This setting captures 
the problem of packet aggregation in convergecasting trees. For this model, a 
2-approximation is known for the variant with deadlines [3], and an O(log^)- 
competitive online algorithm is known for linear delay costs [5], where 9 is the 
the sum of all nodes' ordering costs. 

2 Iterative-Rounding Approximation Algorithm 

In this section we present our Algorithm J MM for JRP-D. The algorithm solves 
the linear program for the relaxation of JRP-D and then rounds the obtained 
fractional solution x to an integral solution. The rounding process is based on 
the method of randomized iterative rounding. We show that Algorithm JMM 
computes an order schedule x whose expected total cost is bounded by e/(e — 1) 
times the optimal cost. 

LP relELxation. Without loss of generality, the orders are placed only at starting 
points of some demand periods. Let D be the given set of demands, and let 
T = {r : (p, r, d) G D} be the set of the starting points for all the demands 
in the instance. The following linear program is the fractional relaxation of the 
natural integer program for JRP-D. 

minimize cost{x) ^ EteT (^^« + Epe{i,...m} ^p^^H 

subject to Xt > x1 for alH e r, p e {1, . . . m} (1) 

T.teT,r<t<d < > 1 for all (p, r,d)eD (2) 

Xt, x^ >0 for alH e r, p e {1, . . . m} 

Let vector x be an optimal fractional solution to this LP relaxation. We show 
how to construct a feasible integral vector x such that cost{x) < ■ cost(x). 

Real time and virtual time. Let us start with some intuition. We will 
use fractional orders as an alternative notion of time. Any real time t € T 
will be associated with the virtual warehouse time interval between X^ter t<T -^t 
and J2t£T t<T -'^t- The concept of virtual time of a retailer p can be defined 
analogously. Since warehouse orders are placed more often than retailers' orders 
(as encoded in constraint (1)), this virtual time will flow at least as fast at the 
warehouse as it flows at any retailer. To guarantee feasibility of the obtained 
solution, it suffices that orders of each retailer happen sufficiently often in his 
virtual time measure (as encoded in constraint (2)). 

Let V = J2tGT -^t denote the (possibly fractional) number of warehouse orders 
in the solution x. Without loss of generality we can assume that the instance 
contains at least one demand, so w > 1. For a € [0, u] define 



rt{a) = min jr e T : Y.teT,t<T > a 



which is a function that translates the warehouse virtual time into real time. For 
a G [0,v] and p G {!,■•■ in} define the quantity 



Density p{a) = ^ , 

which relates different virtual time notions. If Xrt{a) = 0, then x^^j.^^ = as well, 
and we set Densityp{a) = in such case. By constraint (f), Densityp{a) < 1 
for all p and a. Note that in these terms, the (possibly fractional) number of 
orders of retailer p in the solution x is equal to X^tgT ~ Jo ^s'^'S*^2/p(2/) dy- 
We denote this amount by v^. 

Algorithm JMM. In our construction, we will use a random number genera- 
tor Rand{), which generates random numbers from (1/e, 1) with a distribution 
defined by the density function 1/y for y G (1/e, 1). Thus the probability that 
RandQ G {z,z'), for 1/e < z < z' < 1, is equal ln(z'/2;). The expected value of 
this distribution is E[Rand{)] = (e — l)/e. 

To compute x, Algorithm JMM first computes the warehouse ordering times, 
and then, for each retailer p, it determines which orders should p join. The details 
of the algorithm are given in Pseudocode 1. In the warehouse virtual time, its 
ordering times are computed (Lines 2-5) as a sequence {aj)j, with the length 
of each interval [uj, a^+i] determined by RandQ. Then we consider each retailer 
p separately. For each p, we compute its order times (in the warehouse virtual 
time) as a sequence (/3f )i, which is a sub-sequence of {aj)j- The idea is to spread 
these times as sparsely as possible (to minimize cost), while ensuring feasibility. 
After computing each (3^ , the next order cannot be later than at the warehouse 
virtual time 7 such that in the virtual time of retailer p the interval [/3f,7] has 
length 1, that is /^p Density{y) dj/ = 1. So we choose as /3f the last warehouse 

ordering time before 7 (Line 9-11). 



Pseudocode 1 Algorithm JMM 

1: a; -f- 0, Qo 0, j 0, u Ylt<=T ^* 

2: while < u — 1 do 

3: Qj+i + RandQ) 

4: ^rt(aj+i) ^ 1 

5: j ^ J + 1 

6: for p = 1, . . . , m do 

7: /?f'? ^ 0, j 

8: while J^p Density ^(y) Ay > 1 do 

9: 7 <~ min{2 : J^p Densityp{y) dy = 1} 

10: ^ max{Qj : < 7} 

12: i^i + 1 

13: return x 
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Fig. 2. An illustration of virtual times, the algorithm, and the proof of Lemma 4. 
Variables xt are the lengths of the intervals on the warehouse virtual time axis. Vari- 
ables x'^ are represented by the areas of shaded rectangles. Function Density p{) is then 
represented by the upper envelope of these rectangles. 



Let q be the number of warehouse orders generated. So their virtual times 
are ai,a2, ■■■,aq, where v — I < aq < v. It could happen that we will set the 
same Xt to 1 twice (that is, when rt(pij) — rt{aj+i) = t), if Xt > 1/e, but not 
more than twice. This does not affect our analysis. The same comment applies 
to each retailer orders, of course. 

Lemma 1. Algorithm JMM terminates in polynomial time and it returns a fea- 
sible solution X. 

Proof. Let n = \D\ be the total number of demands. The loop in Lines 2-5 
terminates, because in each iteration it increases the current value of aj by at 
least 1/e. Since v < \T\, this loop makes q < \T\/e < n/e iterations. Consider 
now the loop in Lines 8-12, for a fixed retailer p. Since 1 = J^p Densityp{y) dy < 

Jjp 1 dy = 7 — /3f , we have 7 > /3f -I- 1, so there exists a j for which /3f < aj < 7. 
Hence, by the choice of /3f_^_]^ in Line 10, we have > /3f and the algorithm 
moves forward in each loop iteration. As (/3f)i is a sub-sequence of {ctj)j, the 
number of these iterations is at most n/e. The integral in Line 9 is easy to 
compute incrementally in amortized constant time. So overall, the running time 
is 0{nm). 

Finally, we have /^p^^ Density p{y) dy < JJp Density p{y) dy = 1. This implies 
that retailer p does not have a demand whose period is strictly contained in- 
between rt{P^) and ri(/3f^-^), because this demand could not be fully served by 
the fractional solution x. Therefore, placing orders for retailer p in both rt{(3^) 
and rt{j3'^j^-^) suffices to satisfy any demand of p with period intersecting the 
interval [rt(/3f),rt(/3f+i)]. □ 

Number of warehouse orders. For proving the approximation guarantee of 
JMM, we first relate the number of warehouse orders in the solution returned 



by the algorithm to v, the number of fractional warehouse orders in the optimal 
fractional solution x. 

By the choice of the random generator RandQ, the expected distance between 
two consecutive warehouse orders (measured in the virtual warehouse time) is 
E[aj+i— Oj] = (e— 1) /e, and orders are placed while the condition in Line 2 holds. 
Intuitively, the number of warehouse a 's is then around w-e/(e— 1). To formalize 
this intuition, we use the following generalized version of Wald's identity (see, for 
example, [11]), customized to our need. Recall that, given a sequence {Ri)i>i of 
random variables, an integer- valued random variable A'^ > 1 is called a stopping 
time for {Ri)i>i if, for any k, the event 'W — fc" is independent of Rj for j > k. 

Lemma 2. Let i?i,i?2,... be a sequence of (possibly dependent) random vari- 
ables, where Ri € [0, 1] and E[i?i] > ^ for all i. Let N be a stopping time for this 
sequence with finite expectation. Then, E[^^-^ > ^•E[A'^]. Both inequalities 
can be replaced by equalities. 

Lemma 3. The expected number of warehouse orders generated &?/ JMM is at 
most -v, where v is the number of warehouse orders in the optimal fractional 
solution X. 

Proof. The random variables Rj = aj — ctj-i, for j — 1, ...,q, satisfy conditions 
of Lemma 2 with ^ = 1 — 1/e and with the stopping time N — q defined as the 
minimum index for which aq > v ~ 1. Then, 

(l-l/e).EM=EE,L,i?,] 

= m, + Rj] 
< 1 + (w - 1) = -y. 

We thus obtain E[q] < ^ - v. □ 

Number of retailer orders. To obtain an analogous bound on the number 
of retailer p's orders, we need to estimate the frequency of p's orders in x. 
Analogously, to apply Wald's identity, we need to show that when considering 
retailer p in the algorithm (cf . Line 6 and the following) , the expected differences 

between consecutive /3'"s are large enough, i.e., that E[J^p+^ Densityp{y) dy] > 
(e — l)/e. However, for the reason we pin down in the proof of Lemma 4, this 
relation might not be true for the last value of f^f. To aleviate this problem, we 
make the following thought experiment: we consider a modified version of the 
algorithm in which the condition aj < v—l'm Line 2 is replaced by aj < v. This 
may incur at most two additional loop iterations in Lines 2-5 (calls of Randi) 
producing at most two more warehouse orders) . As is a deterministic function 
of the sequence a, the newly produced sequence of (denoted P^) is identical 
to the original one, except that the last retailer order might be issued a little 
later. Thus for analyzing the (expected) number of such orders, we may analyze 
sequence instead of /?''. 



Lemma 4. For any retailer p and two consecutive values I3l and /3f+i, it holds 
that 



E 



+1 

Density [y] dy 



>(e~l)/e. 



Proof. Consider the value of 7 computed in Line 9 in the same iteration which 
computes the value of /3f^_]^. Since 

M+i n 

Density Jy) dy ^ 1 - Density Jy) dy > 1 - - |3^_^_J] , (3) 

it is sufficient to show that £[7 — /3,f_^i] < 1/e. Define p to be the smallest index 
for which Up > j — 1 (cf. Figure 2). 

If 7 < Qfp + 1/e then /3f^]^ = ap with probability 1 (by the definition of 
RandQ), so £[7 — /3f_|_]^] < 1/e is trivial. Thus, in the rest of the proof we can 
assume that 7 > ap + 1/e. 

We now focus on the distance /i = 7 — ofp and consider the random choice of 
the next warehouse order moment ckp+i, which was computed by the algorithm 
as Qfp+i ^ ap + y, where y is the random variable of this particular execution of 
Rand(). (Note that our modification of the sequence a gurantees that ap is not 
the last order, and hence ap+i is well-defined.) 

By the case assumption we have 1/e < /i < 1. In case ckp+i > 7, which is 
equivalent to y > fi, the next retailer p's order is placed in 13^^^ = ap. Then 
7 — f3^_^i = H- In the remaining case, we have a^+i < 7, which is equivalent to 
y<H. Then ^'^^^ > a^+i, so 7-/3f^i < 7-ap+i = 7 - (ap + y) = /x-y. Hence: 



dy 



E[7-,9f+i]< r {^i-y)y-'dy+ [ fiy~' 

Jl/c J II 

= p. [ y~^dy~[ ldy = 1/e, 

J 1/0 Jl/c 

thus completing the proof. □ 

Lemma 5. For any retailer p, the expected number 0/ its orders generated by 
JMM is at most -v^ , where is the number of retailer orders in the optimal 
fractional solution x. 

Proof. For the sequence /3f , , • • • , /3f generated by the algorithm, it holds that 

/ Density Jy) dy ^ / Density Jy) dy - / Density Ay) dy < - I , 
Jq Jo Ji}^^_^ 

where the inequality follows by the while condition in Line 8. Consider now 
the modified sequence — P^, fSl^, . . . , (3^ . Recall that /3f — for i < s 
and 0P > We may now apply Wald's identity: the random variables Rj = 

Jpp Density (y) dy for j ~ 1, s satisfy conditions of Lemma 2 with ^ > 1 — 1/e 



and with the stopping time N — s defined as the minimum index for which 
/g ° Density pljj) dy > w'' — 1. Then, 

(l-l/e).E[s] <EE;=ii?,] 

= E[i?, + E;i!i?,] 

< 1 + K - 1) = yP, 

We thus obtain Efsl < ■ . 

Approximation guarantee. By Lemma 1, Algorithm J MM returns a feasible 
solution X. Furthermore, by Lemma 3 and Lemma 5 the expected number of 
orders (respectively, warehouse and retailer ones) it produces is at most 
times the number of orders in the optimal fractional solution x. Hence, we obtain 
the following bound. 

Theorem 1. Algorithm JMM is a polynomial-time randomized -^^-approxima- 
tion algorithm for JRP-D. 

3 APX-Hardness 

Let JRP-De4 be the restriction of JRP-D where each retailer has at most four de- 
mands and all demand intervals are of the same length. We show that JRP-De4 
is APX-hard, using the result by Alimonti and Kann [1] that Vertex Cover is 
APX-hard even for cubic graphs. Roughly speaking, on the basis of any cubic 
graph G = (V, E) with n vertices (where n is even) and m = l.bn edges, in poly- 
nomial time we construct an instance Jg of JRP-De4, such that the existence 
of a vertex cover for G of size at most K is equivalent to the existence of an 
order schedule for Jg of cost at most 10. 5n -\- K -\- Q (and such equivalence can 
be shown by polynomial-time mappings between solutions). 

Construction of instance Jg- In the following, we fix a cubic graph G 
with vertices vq, . . . ,Vn-i and edges eg, . . . , em_i. We construct Jg consisting 
of 1 -I- m -f 3n retailers, where all the ordering costs are 1, that is L = 1 and 
bp = 1 for all p's. Each retailer has at most four demands and all demand periods 
are of length 4m. 

With each edge Cj we will associate a pair of time points 2j,2j + 1, each 
corresponding to an endpoint of Cj. For each vertex Vi, let e^g^j), e^-j^ji) and 
ej2(i) be the three edges incident to Vi. Vertex Vi will be associated with four 
time points. One of these points is j3i = 8m — i, and the other three points are 
those that represent Vi in the point pairs corresponding to edges. Specifically, 
for a £ {0, 1, 2}, if ej^(i) = (ui, w^') then let ai^a = '^ja{i) if * < otherwise let 
Oii,a = '^jaii) + 1- These three other points are then ai^, a^^i and ai^2- 

Jg will consist of 1 + n -t- m gadgets: a vertex gadget YGi for each vertex Vi, 
an edge gadget EGj for each edge Cj, and one special support gadget SG. We 
now describe these gadgets (see also Figure 3 for reference). 
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Fig. 3. The construction of instance Jc- The figure shows the support gadget SG, an 
edge gadget EGj and a vertex gadget VGi. Shaded regions represent retailers. Thick 
horizontal intervals represent demand periods. 



The support gadget. SG consists simply of one retailer with three demands 
(0, —4m — 1,-1), (0, 2m, 6m) and (0, 8m + 1, 12m + 1). Retailer does not have 
any other demands. Since the demand periods of retailer are disjoint, they will 
require three orders. We will show later that these orders can be assumed to be 
at times —1, 4m and 8m + 1. 

Edge gadgets. For any edge ej, EGj consists of retailer j + 1 with demands 
(j + 1, 2j + 1 - 4m, 2j + 1) and (j + 1, 2j, 2j + 4m). Retailers 1, 2, m will not 
contain any other demands. The demands of each such retailer j + 1 can be 
satisfied with one order in the interval [2j, 2j + 1]. If there is no such order, two 
orders are needed. 

Vertex gadgets. For any vertex v^, VGi will involve three retailers pi^a = 
m + 1 + 3i + a, for a = 0, 1, 2, and their demands. Each retailer pi^a will have 
four demands: 

{Pi,a, Cti^a - 4m, ai^a), {Pi,a, di^aja^^a + 4m), 

{Pi,a,f3i - 'lm,/3i), and {pi^a, Pz, f3i + 4:m). 

For ease of reference, we will denote the periods of these four demands as Qi a, 
Qi,ai Qi,a ^nd Qf^, in the order listed above. Note that n Qi,a = {"i,a}: 
Qi,a ^ Qla = {P^ " 4^, a,^a + 4m] ^ Qf^ O = and that other pairs of 
demand periods of retailer ^ are disjoint. The important property of retailer 
pi^a is that it requires two orders, and that the only way to satisfy his demands 
by two orders is when these orders are at ai^a and Pi. There are infinitely many 
ways to satisfy pi^aS demands with three orders. In particular, orders at times 



— 1, 4m and 8m + 1 work, and these times will be used in our construction for 
this purpose. 

Lemma 6. IfG has a vertex cover U of size K, then there is an order schedule 
for Jg of cost at most 10. 5n + K + 6. 

Proof. To construct a schedule for J'g, we first include warehouse orders at times 

— 1, 4m and 8m + 1, and each of them is joined by retailer (of gadget SG). 
This incurs the cost of 6. 

Next, consider a vertex Vi and the associated gadget VGj. li Vi ^ U, have 
retailers ptfl, pi^i and pi^2 join warehouse orders at times —1, 4m and 8m+l. The 
additional cost for VGj is then 9. If Vi G U, for each i = 0, 1, 2 we make two orders 
from retailer pi^a, one at a^.a and one at f3i. Since the warehouse order at f3i is 
shared between these three retailers, the additional cost for VG, is 1 + 3(1 + 2) = 
10. As is a vertex cover, the choices we made for vertex gadgets imply that, 
for each edge ej, there is now a warehouse order at a time t E {2j, 2j + 1}, that 
corresponds to a vertex that covers ej. We can thus have retailer j + 1 of EGj join 
this order at cost 1. Thus, the total cost is 6+9{n—K) + 10K+m = 10.5n+K+6, 
as required. □ 

Recovering a vertex cover from an order schedule. We now want to show 

that from any order schedule of cost 10. 5n + K + 6 for J'g we can compute a 
vertex cover of G of size K. This part takes some additional work, as we would 
like to assume that the given order schedule S is in the following normal form: 

(nfl) S has warehouse orders at times —1, 4m and 8m + 1, and these orders are 

used by retailer 0. 

(nf2) For each edge gadget EGj, S has an order at a time t G {2j, 2j + 1} that 

involves retailer j + 
(nf3) For each vertex gadget VGi, one of the following two conditions holds: 

— S has orders of retailer pi^a at times /3j and ai^a, for each a = 0, 1, 2. 

— S has orders from retailer pi^a at times —1, 4m and 8m + 1, and there 
are no warehouse orders at times /3i and ai,a) for each a = 0, 1, 2. 

(nf4) There arc no other orders in S. 

The following lemma states that we may assume an order schedule to be in 
the normal form, without loss of generality. 

Lemma 7. For any order schedule S for Jg, we can compute a normal-form 
schedule S' whose cost is at most the cost of S. 

Proof. We prove this claim by gradually modifying S to obtain 5'. We can 

assume that orders occur only at endpoints of demand intervals. We can also 
assume that the first order is at the earliest deadline and the last order is at the 
latest release time. Therefore we can assume that 5 has orders at times —1 and 
8m + 1. Because of demand (0,2m, 6m,), there must be an order between 2m 
and 6m, let's say at time t. There are no demand intervals ending in [2m, 4m), 
so if 2m <t< 4m, we can shift the order at t rightward to 4m. Similarly, there 



are no demand intervals starting in (Am, 6m], so if 4m < t < 6m then we can 
shift the order at t leftward to 4m. This does not increase the cost. This way we 
will have an order at time 4m and no other orders in the interval [2m, 6m]. The 
resulting schedule satisfies Property (nfl). 

Next, consider any edge gadget EGj. If retailer j + 1 is not involved in an 
order at time 2j or 2j' + 1, then it must be involved in two orders, one before 
2j and the other after 2j + 1. We can remove him from these two orders and 
place an order at time 2j (or 2j + 1) without increasing the cost. The resulting 
schedule satisfies Property (nf2). 

Finally, consider any vertex gadget VGi. Suppose that S has a warehouse 
order at time Pi. Let a S {0, 1, 2}. We can assume that retailer pi^a participates 
in the order at time Pi, since otherwise we can replace his order at any later 
time (there must be one in the demand period by the order at Pi. Suppose 
that S does not have an order from retailer pi a at time ai^a- Then retailer pi^a 
must participate in at least two orders in addition to that at Pi . We can remove 
him from these two orders and add an order from pi ^ at time ai^a (or have pi^a 
join the warehouse order at time a^.a if it already exists). This operation does 
not increase the cost. We have thus shown that if there is a warehouse order at 
time Pi then, without loss of generality, we can assume that each retailer p^^a, 
for a ~ 2, has orders at times ai^a a-nd Pi, and no other orders. 

Suppose that S does not have an order at time Pi, and that for some a G 
{0, 1, 2} it has an order at time a^.a that does not involve pi^a- In this case pi^a 
must be in three orders. We can remove him from these orders, add him to the 
warehouse order at a^^aj a-nd add a new order Pi. The cost of the added orders 
is at most 3, so the overall cost will not increase. By the previous paragraph, we 
can then assume that for each of the other two retailers pi.a' of VG^, a' ^ a, he 
also has exactly orders, at ai^a' and Pi. 

The two paragraphs above show that either we can modify S so that the 
orders from VGi satisfy the first condition in Property (nf3), or there are no 
warehouse orders at times Pi, ai^o, a^^i and ai_2- In the latter case though, each 
retailer pi^, pi^i, pi 2 must be in three orders, which, without loss of generality, 
are those at times —1, 4m and 8m + 1, i.e., VG^ satisfies the second condition 
in Property (nf3). 

As the orders described in the proof above are the only necessary orders. 
Property (nf4) holds as well and the lemma follows. □ 

Lemma 8. On the basis of an order schedule S for Jq of cost lO.bn + K + 6, 
we can construct in polynomial time a vertex cover of G of size K. 

Proof. By Lemma 7, we can assume that S is in the normal form. For such 
schedules it is easy to compute the cost. Suppose that S has i orders in the 
interval (6m, 8m], that is at times Pi. For each Vi for which S has an order at 
Pi, we pay 1 for the warehouse cost at Pi, for each edge Cj^t^i) incident to Vi 
retailer pi^a pays cost of 2 and we pay 1 more for the warehouse cost at ai^a- So 
the total cost associated with such w^'s is 10€. For each i for which S does not 
have an order at Pi we pay 9 for the retailer cost (3 for each pi^a)- So the cost 



associated with such u^'s is 9{n — £). We then have additional cost m = 1.5n for 
the retailer cost of the gadgets EG^ , plus 3 for the warehouse cost at —1, 4m, 
8to + 1 and 3 for the retailer cost for gadget SG at these times. This gives us 
total cost 10.571 + £ + 6, and we can conclude that we must have £ — K. 

Define now U to be the set of those vertices Vi for which S makes an order 
at time /3i from the retailers in VG^. Consider some edge Cj. Using the normal 
form, we first obtain that there is an order at time 2j or 2j + l, say at 2j. Ghoose 
i and a for which 2j — ai,a- By the normal form, gadget VG^ satisfies the first 
condition in Property (nf3), which in turn implies that Vi € U. This shows that 
each edge ej is covered by U, so G has a vertex cover of size K. □ 

Theorem 2. Problem JRP-De^, that is the restriction of JRP-D to instances 
with equal demand intervals and at most four demands per retailer, is AVlL-hard. 

Proof. By [1], there exists an e > 0, such that there is no polynomia-time (1 + e)- 
approximation algorithm for the vertex cover problem in cubic graphs, unless 
P = NP. Towards a contradiction, assume that we can compute in polynomial- 
time a (1 + e/24)-approximation for JRP-De4. Fix any cubic graph G of n > 6 
vertices and m edges with a vertex cover of size K and construct an instance 
J^G- By Lemma 6, J7g has an order schedule of cost at most 10.5n + K + 6, so our 
approximation algorithm finds a schedule of cost C < (l-|-e/24) • (10.5n-|-i^-|-6). 
As the graph is cubic, K > m/3 = n/2, and therefore 10.571 + 6 < 11.571 < 23i^. 
Thus, C < 10.5n + {l + e)K+6. (In fact, as all the costs are integers, the obtained 
order schedule has cost at most 10.571 + [(1 + e)K\ + 6.) Finally, we may use 
the construction of Lemma 8 to compute a vertex cover for G of size at most 
[(1 + e)K\ , which contradicts the result of [1]. □ 

4 1.5- Approximation for Equal-Length Demand Periods 

In this section, we present a 1.5-competitive algorithm for the case where all the 
demand periods are of equal length. Without loss of generality, we will assume 
that the length of each period is 1. 

We denote the input instance by I. Let the width of an instance be the 
difference between the deadline of the last demand and the release time of the 
first one. The building block of our approach is an algorithm that creates an 
optimal solution to an instance of width at most 3. Later, we divide I into 
overlapping sub-instances of width 3, solve each of them optimally, and finally 
show that by aggregating their solutions we obtain a 1.5-approximation for I. 

Lemma 9. A solution to an instance J of width at most 3 consisting of unit- 
length demand periods can be computed in polynomial time. 

Proof. We shift all demands in time, so that is entirely contained in inter- 
val [0, 3]. Recall that L is the warehouse ordering cost and bp is the ordering cost 
of retailer p € {l,2,...,m}. Without loss of generality, we can assume that all 
retailers 1, ..,771 have some demands. 



Let dmin be the first deadline of a demand from J and rmax the last release 
time. If Tujax < c'minj then placing one order at any time from [rmaxi c^min] is 
sufficient (and necessary). Its cost is then equal to i + hp. 

Thus, in the following we focus on the case dmin < ^max- Any feasible solution 
has to plac;c an order at or before dmin and at or after rmax- Furthermore, by 
shifting these orders we may assume that the first and last orders occur exactly 
at times dmin and Tmax, respectively. 

The problem is thus to choose a set T of warehouse ordering times that 
contains dmin, "Tmaxj and possibly other times from the interval (dmin, ''max), and 
then to decide, for each retailer p, which warehouse orders it joins. Note that 
''max — dmin < 1, and therefore each demand period contains dmin, ?'max, or both. 
Hence, all demands of a retailer p can be satisfied by joining the warehouse 
orders at times dmin and rmax at additional cost of 26p. It is possible to reduce 
the retailer ordering cost to hp if (and only if) there is a warehouse order that 
occurs within Dp, where Dp is the intersection of all demand periods of retailer 
p. (To this end. Dp has to be non-empty.) 

Hence, the optimal cost for J can be expressed as the sum of four parts: 

(i) the unavoidable ordering cost hp for each retailer p, 

(ii) the additional ordering cost hp for each retailer p with empty Dp, 

(iii) the total warehouse ordering cost L- |T|, and 

(iv) the additional ordering cost hp for each retailer p whose Dp is non-empty 
and does not contain any ordering time from T. 

As the first two parts of the cost are independent of T, we focus on minimizing 

the sum of parts (iii) and (iv) that we call the adjusted cost. Let C{x) be the 
minimum possible adjusted cost "in the interval [dmim x]" under the assumption 
that there is an order at time x. Formally, C{x) is the minimum, over all choices 
of sets T C [dmin,a^] that contain dmin and x, oi L ■ \T\ + X]peQ(r) ^P' "where 
Q{T) is the set of retailers p for which Dp^% and Dp C [0,a;] - T. (Note that 
the second term consists of expenditures that actually occur outside the interval 

[dmin,-^;]-) 

As there are no Dp's strictly to the left of dmin, C'(dmin) = L. Furthermore, 
to compute C{x) for any x G (dmin, ''max], we can express it recursively using 
the value of C{y) for y e [dmin,-^) being the warehouse order time immediately 
preceding x in the set T that realizes C{x). This gives us the formula 

C{x)=L+ rnin (c(y) + ^ bp) . 

In the minimum above, we may restrict computation of C{x) to .t's and j/'s 
that are ends of demand periods. Hence, the actual values of function C(-) can 
be computed by dynamic programming in polynomial time. Finally, the total 
adjusted cost is equal to C(rmax)- Once we computed the minimum adjusted cost, 
recovering the actual orders can be performed by a straightforward extension of 
the dynamic programming presented above. □ 



Now, we construct an approximate solution for the original instance X con- 
sisting of unit-length demand periods. For i e N, let I,j be the sub-instance 
containing all demands entirely contained in [i, j -I- 3). By Lemma 9, an optimal 
solution for X^, denoted can be computed in polynomial time. Let Sq be 

the solution created by aggregating A(Xq) , A(Xi) , A{^^), . . . and S\ by aggregat- 
ing ^(Ti), ^(Ta), ^(Xs), . . .. Among solutions 5*0 and 5*1, we output the one with 
the smaller cost. 

Theorem 3. The above algorithm produces a feasible solution of cost at most 
1.5 times the optimum cost. 

Proof. Each unit-length demand of instance X is entirely contained in some l2fe 
for some fc € N. Hence, it is satisfied in A{X2k), and thus also in So, which yields 
the feasibility of 5*0. An analogous argument shows the feasibility of 5*1. 

To estimate the approximation ratio, we fix an optimal solution Opt for 
instance I and let opt^ be the cost of Opt's orders in the interval [i, i + 1). Note 
that Opt's orders in [i, i + 3) satisfy all demands contained entirely in [i, i + 3). 
Since A{Xi) is an optimal solution for these demands, we have cost{A{Ii)) < 
opti + opti+i -f- opti+2 and, by taking the sum, we obtain cost{So) + cost{Si) < 

cost{A{Ii)) < 3 • cost(OPT). Therefore, either of the two solutions (5*0 or Si) 
has cost at most 1.5 • cos<(Opt). □ 

5 Final Comments 

The main result we presented is an e/(e — l)-approximation algorithm for JRP-D. 
The main open question is of course whether and by how much the approxima- 
tion ratio for JRP-D can be improved, and what ratio can be achieved with a 
deterministic algorithm. A careful reader might have noticed that some improve- 
ments should be still possible: since the density function of the next order used 
in our algorithm is non-zero starting at 1/e < 0.5 (in terms of the fractional 
solution "mass"), for time points that are farther than 2/e and not farther than 
1 from the last warehouse order, with some small probability our algorithm gets 
a second chance to use them for the next order. We have not pursued this ap- 
proach, since the improvement in the ratio will be very minor and its analysis 
looks tedious. 

There is a simple algorithm for JRP-D that provides a (1, 2)-approxiniation, 
in the following sense: its warehouse order cost is not larger than that in the 
optimum, while its retailer order cost is at most twice that in the optimum [10]. 
One can then try to balance the two approaches by choosing each algorithm with 
a certain probability. This simple approach does not improve the ratio. But it 
may be possible to achieve a better ratio if, instead of using our algorithm as 
presented, we appropriately adjust the probability distribution. 

If we parametrize JRP-D by the maximum number p of demand periods of 
each retailer, its complexity status is essentially resolved: for p > 3 the problem 
is APX-hard [10], while for p < 2 it can be solved in polynomial time (by a greedy 
algorithm for p = 1 and dynamic programming for p = 2) . In case of equal-length 



demand periods, we showed that the problem is APX-hard for p > 4, but the 
case p — 3 remains open, and it would be nice to settle this case as well. We 
conjecture that in this case the problem is NP-complete. 
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